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Abstract 

In this paper, we consider a Fisher-KPP equation with an advection 
term and two free boundaries, which models the behavior of an invasive 
species in one dimension space. When spreading happens (that is, the 
solution converges to a positive constant), we use phase plane analysis and 
upper/lower solutions to prove that the rightward and leftward asymptotic 
spreading speeds exist, both are positive constants. Moreover, one of 
them is bigger and the other is smaller than the spreading speed in the 
corresponding problem without advection term. 



1 Introduction 

In 2010, Du and Lin [5] studied the following Fisher-KPP problem with free 
boundaries: 



u t — du xx = u(l — u), g(t) < x < h(t), t > 0, 

u(t,g(t))=0, g'(t) = -tiu x (t,g(t)), t>0, 

u(t,h(t))=0, h'(t) = -fiu x (t,h(t)), t>0, 

—g(0) — h(0) — Hq, u(0,x) — uq(x), —ho < x < ho, 



(Po) 



where d and /i are positive constants, the initial function uo(x) satisfies 

u e C 2 ({—h ,h }), u (±h ) — and uq > in (— ho, ho), (1.1) 
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for some ho > 0. They used QPq| ) to model the spreading of a new or invasive 
species with population density u(t, x) over a one dimensional habitat, with the 
free boundaries x = g(t), h(t) representing the expanding fronts. They obtained 
a dichotomy result, that is, either spreading happens (u(t, •) — > 1 locally uni- 
formly in K and h(t), —g(t) — > oo as t — > oo) or vanishing happens (u(t, ■) — > 
uniformly in [g(t),h(t)] as i — > oo and /i(t) — g(i) < oo). Furthermore, when 
spreading happens, they obtained the existence of the asymptotic spreading 
speed (0 Proposition 4.1]): 

c * := Um 'M = i im ZiQl > o. (1.2) 

t— voo t t— >oo t 

Recently, further extensions have been done, for example, Du and Guo [H [S] 
studied the problem in higher dimension spaces and in heterogeneous environ- 
ment. Du and Lou [7] studied the problem with general nonlinear /, including 
general monostable, bistable and combustion types of /. Among others, they 
all proved that the asymptotic spreading speed when spreading happens is the 
same positive constant in any direction. 

However, some species prefers to move towards one direction because of 
rich resource, appropriate climate, etc. Some diseases spread along the wind 
direction. In 2009, Maidana and Yang in [9] studied the propagation of West 
Nile Virus from New York City to California state. It was observed that West 
Nile Virus appeared for the first time in New York City in the summer of 1999. 
In the second year the wave front travels 187km to the north and 1100km to 
the south. Therefore, they took account of the advection movement and showed 
that bird advection becomes an important factor for lower mosquito biting rates. 
Recently, Averill in [1] considered the effect of intermediate advection on the 
dynamics of two-species competition system, and provides a concrete range of 
advection strength for the coexistence of two competing species. Moreover, 
three different kinds of transitions from small advection to large advection were 
illustrates theoretically and numerically. 

What is the difference between the asymptotic spreading speed of the left 
frontier and that of the right frontier when invasive species is spreading? To 
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address the question, in this paper we study the following problem with an 
advection term: 

u t - du xx + f3u x = u(l - u), g{t) < x < h(t), t > 0, 

u(t,g{t)) = 0, g'(t) = -fiu x (t,g(t)), t > 0, 

u(t, h(t)) = 0, h'(t) = -fm x (t, h(t)), t > 0, K l > 

— 5(0) = MO) = ^0: u(0,x) = Uo(x), —hg <x< ho, 

where d, fx, ho, uo are as above and /3 > is a constant. 

By a similar argument as in [HO [7], we have the following basic results. 

(i) Problem has a time global solution (u, g, h) with u G c*i+" A2+a(( 0; 00) 
x /i(t)]) and g, h G C 1+Q / 2 ([0, 00)) for any a G (0, 1); 

(ii) < u(t,x) < Ci for g(t) < x < h(t), t > and < -g'(t),h'(t) < C 2 for 
t > 0, where C\ and C2 are constants independent of t. 

In a forthcoming paper [5] , we studied the asymptotic behavior of the solutions 
of Q_Pi| ). More precisely, we gave some sufficient conditions for spreading and 
some sufficient conditions for vanishing. It turns out that spreading happens 
only if 

< P < 2Vrf. (1.3) 

This paper is devoted to the difference between the leftward and rightward 
asymptotic spreading speeds induced by the advection term (3u x . This is an 
interesting problem from ecological point of view. 

Therorem 1.1 Assume < /? < 2\[d. Let (u,g,h) be a solution of 
( |PiD for which spreading happens. Then the leftward and rightward asymptotic 
spreading speeds exist: 

* ,. -<?(*) * r Kt) 

c ; := hm , c r := lim . 

t— >oo t t— >oo t 

Moreover, < c ; * < c* < c*, where c* is the spreading speed of the solution of 



c* is given in (11.2[) which is nothing but ko in [6j Proposition 4.1], or c* in [7, 
Theorem 1.10]. It depends on the parameter fi. Similarly c[ and c* depend 
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on [i. Clearly, c* and c\ also depend on j3. On these dependence we have the 
following results. 

Therorem 1.2 (i) If ft £ (0,2-y/d) is fixed, then c*,c*,c* are strictly 
increasing in [i, and 

lim c,* = 0, lim c\ = 2Vd - /?, 

fj,— n— >oo 

lim c* = 0, lim c* = 2\Q, 
lim c* = 0, lim c* = 2 + /3; 

(ii) i/ /x zs fixed, then c*, — c ; * are strictly increasing in j3, and 
lim cf = lim c* = c*. lim cf = 0. 

2 Semi-waves and spreading speeds 

Throughout this section we assume that (|1.3[) holds and that (w, g, h) is a solu- 
tion of flPiD for which spreading happens, that is, h(t), —g(t) — > oo (i —> oo), 
and ■) — > 1 locally uniformly in R. Denote /(u) := u(l — u) for convenience. 
We remark that the approaches below remain valid for general monostable non- 
linear /. 

To determine the spreading speed, we will construct upper and lower solu- 
tions based on semi-waves. 

2.1 Phase plane analysis and semi- waves 

We call q(z) a semi-wave with speed c if (c, q(z)) satisfies 



q" - ^q' + ^ = forze(0,oo), 
g(0) = 0, q(oo) = 1, q(z) > for z £ (0, oo). 

The first equation in this problem is equivalent to the following system: 



(2.1) 



„/_c-j9 /(g) ( 2 - 2 ) 
r ~ d V d ■ 

A solution (q(z),p(z)) of this system traces out a trajectory in the q,p-plane or, 
as it is usually called, the phase plane (cf. [H [3J [7J [TU] ) . Such a trajectory has 
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slope 

dp = c-P _ /(g) ^ ^ 

dq d dp 

at any point where p ^ 0. Here we are only interested in a trajectory of (12.21) 
that starts from the point (0, oS) with some uj > and ends at the point (1,0) 
as z — » +oo. 

For any fixed c > 0, (0,0) and (1,0) are critical points of the system (|2.2[) . 
The eigenvalues of the corresponding linearizations are 



, _ c-p± ^{c-Pf-Ad 



A ° = 2d (at ( °' 0)) ' 



c 



/3 ± ^(c - /3) 2 + Ad 



K = W \ d ' (at (1,0)), 

respectively. Thus (1,0) is a saddle point and (0,0) is 

(i) a center or a spiral point, if < c < P + 2Va; 

(ii) a nodal point, if c > /3 + 2vd. 

Therefore, by the theory of ODE (cf. [10]), there exactly two trajectories of 
(|2.2jl that approach (1,0) from q < 1. One of them, denoted by T^, has slope 
< at (1, 0). Suppose that T£ is expressed by a function p — Pfiq)- Then 
p = Pfiq) satisfies (|2.3p and T° lies in the semistrip 

S = {(g,p):0<g<l,p>0}. 

T r c is a trajectory through (1,0) and (0,P r c (0+)) for some P r c (0+) > 0. The 
following are well known results (cf . [3J [JJ [IT] ) . 

Proposition 2.1 Lef c° := 2\Q + p. Then 



(i) for any c £ [0, c°), P^(0) is positive, continuous, strictly decreasing in c, 
and lim P r c (0) = 0; 

(ii) for any c > c% P r c (0+) = 0. 
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In case (ii), each T° is a trajectory in S through (0,0) and (1,0), and so 
it corresponds to a traveling wave with speed c, is nothing but the minimal 
speed of these traveling waves. 

Denote ((c) := Pf(0) - f for c G [0,c°). In view of Proposition CD] <(c) is 
continuous and strictly decreasing in c £ [0,c°), and it satisfies 



C(0) = P°(0) > P/(0) = sl- d f Q f(s)ds > 0, 

c((c?n = ~<o. 

A* 

Thus there exists a unique c; G (0,c°) such that C(c*) = 0, i.e. P r Cr (0) = c*//i. 
Summarizing the above results we have the following proposition. 

Proposition 2.2 Problem (|2.1[) /ias exactly one solution (c,q) = (c*,q*) 
such that 

M<7 r *)'(0) = c;. (2.4) 

Moreover, c* G (0,/3 + 2\/d). 

Later we will use this semi-wave to estimate the rightward spreading speed. 
Similarly, to estimate the leftward spreading speed, we will need another semi- 
wave traveling to left, which is a solution of the following problem: 

q »_c+l q , + m = o forze(0,oo), (2) 
q{0) = 0, q{oo) = 1, q{z) > for z £ (0, oo). V ' ' 



Similar as above, this problem can be studied by considering the problem 

dp_c + fi f(q) 



(2.6) 



dq d dp 

in the g,p-phase plane, where p — q' . Denote Pf(g) the solution of this equation 
whose trajectory through (1,0) and (0, P ; c (0 + )) for some P ; c (0 + ) > 0. In a 
similar way as above we have the following results. 

Proposition 2.3 Let c? := 2y/d- /3. Then 

(i) {or any c G [0, c°), Pf(0) is positive, continuous, strictly decreasing in c, 
and lim if(0) = 0; 



G 



(ii) for any c > c%, P, c (0+) = 0. 

Proposition 2.4 Problem (|2.5I) has exactly one solution (c,q) = (c*,q*) 
such that 

M(<z;)'(0) = c?. (2.7) 

Moreover, c z * G (0, c°). 

Next, we make suitable perturbations of f[u) to derive corresponding semi- 
waves that can be used to construct upper and lower solutions of fl-PiD - 
For any small e > 0, set 

fe(u) ■= f(u) - i~ u2 = u ( 1 ~ Y~~e U ) ' 

£(u) := /(«) + = u(l - ^ti) . 

Note that / e (u) is strictly decreasing in e and it has exactly two zeros and 
1 — e. / E (u) is strictly increasing in e and it has exactly two zeros and 1 + e. In 
a similar way as above, we know that problem (|2.1[) with / replaced by f £ (resp. 
f e ) has exactly one solution (c*,q*) with n(q*)'(0) = c* and c* G (0, c°) (resp. 
a solution (c£,g£) with /i(gj:)'(0) = and c* r G (0,c°)), where c° = 2\fd + (3. 
Similarly, problem (12.5[) with / replaced by f e (resp. f e ) has exactly one solution 
(c*,5f) with /j,(qf)'(0) = Cjf and cf G (0,c°) (resp. a solution (cf,gf) with 
(i(qf)'(O) — c*i and c^ 1 G (0, c?)), where c° = 2v^ - /3. 

Proposition 2.5 T/ie following conclusions hold. 

c < c„ < c r , hm c„ = hm c_ = c„ , 

and 

c; < c ; < Cj , hm c z = hm q = c, . 

Proof. We first prove c* < c*. For any c G [0, c°), consider the problem (|2 .3[) 
with / replaced by / E , denote the solution with trajectory through the critical 
point (0, 1 — e) in the phase plane by P^ £ {q). Similar as Proposition 12. 1 1 (i) we 
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have P r c , e (0) > for all c e [0, c°). Moreover, P r %(g) < P r c (g) (g € (0, 1 - e]) by 
/e(?) < /(g) (0 < q < 1 — e)- We now prove 

< P r c £ (0) < P r c (0) for c G [0, e°). (2.8) 

Otherwise, Pj: E (0) = P,? (0), and so the function 77(g) := P^r (g) — P£ e {q) satisfies 

T]'<a{q)r] (0 < q < 1 — e), r;(0) = 0, 

where a(g) := /^[dP^P^g)]- 1 . This implies that 77(g) < (0 < g < 1-e), 
a contradiction. 

Denote C(c) := P r c !£ (0) - ^ Then (J2U) implies that 

C(c)<C(c) force [0,c°). 

Similar as above, both £(e) and £(c) are continuous and strictly decreasing 
functions in [0,c£), and 

C((c°)-)=C((c°)-) = -i 

M 

Therefore c* < c* by their definitions: C(c*) = C( c *) = 0. 

Next we prove lim £ ^ c* = c*. It is sufficient to show that, for any c € [0, c°), 

P r y0) ->■ P r c (0) as e 0. (2.9) 

By the monotonicity of f e , it is easily seen that P^ E (q) is monotonically decreas- 
ing in e, and it is bounded from above by P° (g). Therefore, as e — > 0, P,? e (g) 
converges to some function P(g) in C 1 ([0,1 — <5]) for any < 6 < 1. Clearly, 
77 = P(g) corresponds to a trajectory of (|2.2p that approaches (1, 0) in the phase 
plane with a non-positive slope at (1,0). Consequently, R(q) = P°(g), and so 
(|23| is proved. 

Other conclusions can be proved in a similar way as above. □ 
2.2 Asymptotic spreading speed 

Proof of Theorem \1.1\ We first estimate the rightward asymptotic spreading 
speed. For any small e > we define 

w(t, x) :— cft(c*rt — x), x € [0, c£i], 
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Since (c*,q*(z)) satisfies 

q" - ^q' + = for z G (0, oo), 
q(0) = 0, g(oo) = 1 - e, q'(0) = £ and <j'(z) > (z > 0), 

we have 

{(5(t,x) < 1 - £, tO t - + /Ju^ = f £ (w) for x G [0,c£t], i > 0, 

and 

fi;(i,c£t) =0, cj; = — fj,w x (t,c^.t) for i > 0. 

Since we are considering the spreading case, we have limt_>oo u(t, •) = 1 
locally uniformly in R. In particular, 

u(t, 0) > 1 - e for t > T 

for some T > 0. Thus (w(t, x),c*t) is a lower solution of flPi| ) on {(t,x) \ x £ 
[0,c*t], t > 0} by comparison principle (cf. [3 [7]), and 

3Jt < fe(t + T), w(t, x) < u(t + T, x) in {(f,£c) | x G [0,cpt], i > 0}. 

This implies that 

lirninf M!l >£. (2.10) 

t— >oo £ 

Next we estimate the upper bound of the rightward spreading speed. Con- 
sider the problem 

V'(t) = f(v) (t>0), »7(0) = ||tio||oo + l. 
A simple comparison shows that 

u(t,x)<ri(t)~ (l- e^Y 1 for *G [<?(*), t > 0. 

Hence for any small e > 0, there exists T > such that 

u(t,x) < 1 + | for z G [0, /i(t)], t > f . 

Recall that (c^,$j!(z)) is a solution of problem (|2.ip with / replaced by / £ and 
g£(oo) = 1 + e. Hence there exists x > h(T) large such that 

u(T,x) < 1 + | <£(s-ss) forxe [0,/i(T)]. 
9 



Define 

w(t, x) := cfi. (c*t + x - x) for x € [0, c^t + x] , t > 0. 

Then (w,c*t+x) is an upper solution of (Fj} on {(t, a;) |x€ [0, ft(i+T)], t > 0}, 
and by the comparison principle (cf. [5J [7] ) we have 

ft(t + T) <c%.t + x, u(t + T,x) < w(t, x) for x e [0, + f )] and t > 0. 

This implies that 

limsup-^- < c^. (2.11) 
Since the limits (|2 . 10[) and (|2 . 1 1 1) hold for any small e > 0, we have 

hm = c r 

t— >oo t 

by Proposition ^. 51 The leftward spreading speed 

hm = c, 

i->oo £ 1 

is proved similarly. 

In [5J[7], the authors considered problem QPq| ), that is, problem QPiD without 
advection term (i.e. /3 = 0). Among others, they showed that the asymptotic 
spreading speed is characterized by the following problem 

dq a dp 

Using a similar phase plane analysis as above, the authors in [Tj proved that 
problem (|2.12p has a solution (c, P c (q)) for each c £ [0,2V")- Moreover, they 
proved that P c (0) \ as c / c° := 2%/d; (|2~T2|) has a unique solution 
(c*,P c *(g)) such that ^P c *(0) = c*. This c* is nothing but the rightward 
and leftward spreading speeds (cf. OE]). 

Combining the above phase plane analysis we have the following conclusions: 

1. P/^O) = P c (0) = P r c +^(0) for all c £ [j3,2Vd); 

2. Pf(0) (resp. P c (0),P^(0)) is continuous and strictly decreasing in c £ 
[0,2Vd-P] (resp. c£ [0,2Vd], c£ [0, 2\/d + /?]). 
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Define three new functions 7 r (c), 7(c) and 7;(c) by 

7r (c) := P r c (0) for c € [0, 2^d + (3), 7(c) := P c (0) for c e [0, 2Vd) 

and 

7l(c) :=i?(0) for ce [0,2^3-0). 

Then, in the c, 7-plane their graphes lie in the first quadrant (see Figure 1) and 
these graphes contact the straight line 7 = — at points (c*,^-), (c* , — ) and 
(c*, ^j), respectively. Therefore, c ; * < c* < c*. This completes the proof of 
Theorem O □ 

Proof of Theorem ] The conclusions in Theorem 11.21 follow from a simple 
analysis on the relations among the graphes of 7 r (c), 7(c), 7/(c) and c//x in Figure 
1. □ 

References 

[1] I. E. Averill, The Effect of Intermediate Advection on Two Competing 
Species, Doctor of Philosophy, Ohio State University, Mathematics, 2012. 

[2] D.G. Aronson and H.F. Weinberger, Nonlinear diffusion in population ge- 
netics, conbustion, and nerve pulse propagation, in Partial Differential 
Equations and Related Topics, Lecture Notes in Math. 446, Springer, 
Berlin, 1975, pp. 5-49. 

[3] D.G. Aronson and H.F. Weinberger, Multidimensional nonlinear diffusion 
arising in population genetics, Adv. in Math., 30 (1978), 33-76. 

[4] Y. Du and Z.M. Guo, Spreading-vanishing dichotomy in a diffusive logistic 
model with a free boundary, II, J. Diff. Eqns., 250 (2011), 4336-4366. 

[5] Y. Du and Z.M. Guo, The Stefan problem for the Fisher-KPP equation, J. 
Diff. Eqns., 253 (2012), 996-1035. 

[6] Y. Du and Z.G. Lin, Spreading-vanishing dichtomy in the diffusive logistic 
model with a free boundary, SIAM J. Math. Anal., 42 (2010), 377-405. 



11 



[7] Y. Du and B.D. Lou, Spreading and vanishing in non- 
linear diffusion problems with free boundaries, preprint 
( http: / /turing. une.edu.au/~ydu/rpub. html | . 

[8] H. Gu, Z.G. Lin and B.D. Lou, A reaction- advection- diffusion equation with 
free boundaries, preprint. 

[9] N. A. Maidana, H. Yang, Spatial spreading of West Nile Virus described 
by traveling waves, Journal of Theoretical Biology, 258(2009), 403C417 

[10] I.G. Petrovski, Ordinary Differential Equations, Prentice-Hall, Englewood 
Cliffs, New Jersey, 1966, Dover, New York, 1973. 

[11] A.I. Volpert, V.A. Volpert and V.A. Volpert, Traveling Wave Solutions of 
Parabolic Systems, Amer. Math. Soc., Providence, Rhode Island, 1994. 



12 



